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Abstract 

The transverse-field Ising model with random exchange interactions in finite dimensions are 
investigated by means of a real-space renormalization group method. The scheme yields the exact 
values of the critical point and critical exponent v in one dimension, and some previous results in 
the case of random ferromagnetic interactions are reproduced in two and three dimensions. We 
apply the scheme to spin glasses in transverse fields in two and three dimensions, which have not 
been analyzed very extensively. The phase diagrams and the critical exponent v are obtained, and 
evidence for the existence of an infinite-randomness fixed point in these models is found. 

PACS numbers: 75.50.Lk, 05.30.Rt, 64.60.F-, 05.10.Cc 



I. INTRODUCTION 



Random quantum spin systems have produced fruitful physics beyond our expectations. 
One of the turning points that opened the door to this exciting topic was the introduction 
of random impurities by McCoy and Wu [1] to the two dimensional classical Ising model. 
The way they introduced the random impurities may seem to be artificial, but using the 
transfer matrix method, one finds a natural quantum spin system equivalent to the McCoy- 
Wu model, the random transverse-field Ising spin chain. After the work by McCoy and Wu, 
a detailed analysis of the corresponding quantum system revealed interesting and peculiar 
features of random quantum spin chains [2J, which are characterized by a so-called infinite- 
randomness fixed point 

Critical phenomena in several low- dimensional random quantum spin systems are con- 
sidered to be controlled by infinite-randomness fixed points in the renormalization group 
picture 4|. As a random system with an infinite-randomness fixed point is coarse-grained, 
the width of distributions grows indefinitely for the logarithms of the parameters in the 
renormalized Hamiltonian. In other words, the randomness is infinitely amplified. The fixed 
point represents broad distributions of physical quantities. As a consequence, rare regions 
in the system, which are usually expected not to affect macroscopic properties, influence 
some of the behavior of the system. Average values of physical quantities, in fact, show 
differences from typical ones. In addition, a characteristic time scale is also influenced, and 
the dynamical critical exponent becomes infinity. 

Such behavior was discovered in the random transverse-field Ising spin chain through 
an analytical study [2j. The same method has been applied to higher dimensions js-ll], 



which numerically verified that higher-dimensional random ferromagnet systems also have 



12| 



in 



the same properties. There is, however, another approach with a conflicting result 
higher dimensions, and hence further investigations are necessary to resolve the controversy. 

The quantum spin glass model is a particularly important issue among random quantum 
spin systems, where the properties of critical phenomena have not been cleared. Although 
the quantum spin glass model has been presumed to have an infinite-randomness fixed point 
based on an analogy with the random ferromagnet model jfj] , strong disorder renormalization 
group (J], which is used in most current activities in random quantum spin systems, has not 



been successfully settled the issue so far. On the other hand, other numerical estimates 
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Q have indicated results against the conjecture, but the system size used in these studies 
might be too small. 

In the present paper, we study the nature of critical phenomena in the random ferromag- 
net Ising model and the spin glass in transverse fields in one, two, and three dimensions with 
a method based on the real-space renormalization group. The method has been developed 
for the non-random model [15|, and it is a variant of the conventional block-spin transforma- 
tion, which has been considered to be a crude approximation in general. Nevertheless, our 
method yields accurate estimates of the critical exponent v in non-random systems 151 ] . We 
apply it to models with randomness in this paper. A special emphasis is placed on whether 
the infinite-randomness fixed point exists in these models. 

Indeed, we have found infinite-randomness fixed points in these models. The critical 
exponent v for the correlation length has been calculated, and its values of the random 
ferromagnet Ising model and the spin glass in transverse fields have been found to be very 
close to each other. This result suggests that the spin glass in transverse fields belongs to 
the same universality class as the random ferromagnet in transverse fields. 

In Sec. [Til we introduce the real-space renormalization group scheme for the transverse- 
field Ising chain with randomness after a short review of the model, and we show that some 
previous results are exactly reproduced under our scheme. In particular, the existence of 
the infinite-randomness fixed point is verified. Our method is generalized to higher spatial 
dimensions in Sees. II III and IIVI and the random ferromagnet model in transverse fields and 
spin glass in transverse fields are investigated. Our study is concluded in Sec. [V] 

II. ONE DIMENSION 

A. Transverse-field Ising spin chain 



12, 



Let us recall a few previous results [2|, |16| related to the random transverse-field Ising 
chain 

i=i i=i 

where a" denotes the a-component of the Pauli matrix on site i, and the boundary condition 
is periodic Oi = <Ti + N, where the number of spins is N, and we assume N is even. The cou- 
plings Jj and transverse fields are random variables independently distributed. Without 



loss of generality, we restrict random variables to be positive. This model already has been 
analyzed in detail j^j], and in the next section we try to reproduce previously known results 
with our real-space renormalization-group approach. 

When the average of transverse fields is much larger than that of couplings, the system lies 
in the paramagnetic phase, that is, the expectation value of a z is zero. A phase transition to 
the ferromagnetic phase takes place at some point as we reduce the average values of fields, 
and the system obtains a finite expectation value of a z . The para-ferro transition point has 
been analytically obtained fl6j |. 

N N 

^logj^^iogr,, (2) 

i=l i=l 



which includes the non-random case 



The model has an infinite-randomness fixed point, and critical phenomena controlled 
by the fixed point has a considerable difference from the conventional ones. One of the 
characteristic features of the critical phenomena appears in the behavior of correlation length 
near the fixed point. There are two kinds of critical exponent v about the divergence 
of correlation length, and the two critical exponents take different values in the random 



transverse- field Ising chain, v — 2 [2] and z/ typ = 1 18|. The first one v is for the (average) 
correlation length £, and it correctly takes into account the effects of randomness. The other 
exponent z/ typ is for the typical correlation length £ typ , which describes the typical behavior 
of the system and does not reflect the influence of rare events on macroscopic properties. 
The correlation length £ properly describing the effect of rare events in the random system 
is defined as the largest length L, where the probability that all the spins in the block of 



length L are correlated, exceeds some finite value [2 



n 



In the random transverse-field Ising model, average and variance of (log r — log J) plays 
an important role in the determination of the correlation length exponent Q]. We now 
express the average as A, 

A = -Y^og-. (3) 

i=l 1 

According to Eq. (j2J), the phase transition in one dimension occurs when this average is 
equal to zero. The system lies in the paramagnetic phase if A > , and the system in the 
ferromagnetic phase has a negative A. We have to treat the average per block of length L 



to find the correlation length as the definition explained above. The average for the block is 

' 1=1 ieb/ 

where bi denotes the set of spins in the block labeled J, and N/ L is the number of blocks. 
Roughly speaking, the variance of (log r — log J) for a block is LV, where V expresses the 



variance for a spin. If LA < yLV, all spins are correlated in many blocks. Otherwise there 
are few blocks containing perfectly correlated spins. The relation £A ~ V£V is thus a good 
estimate for the correlation length. Hence, we have 



e~(-4) • (5) 



The correlation length diverges around the transition point A = 0, and we find A/yV to be 
a proper parameter measuring the "distance" from the critical point. In general, including 
the case A c ^ 0, we ought to expect the relation between the correlation length £ and the 
critical exponent v to be represented as 



The above-mentioned crude estimate suggests v — 2. If we measure the distance without 
variance, or the effect of rare events, we can associate the typical correlation length £ typ with 
the critical exponent z/ typ , 

t, tyv ~{A-A c y v ^. (7) 



B. Real-space renormalization group in one dimension 

We develop a real-space renormalization group procedure to the random transverse-field 
Isine model in one dimension Eq. ([IT) at zero temperature. This is a natural generalization 
of our previous method [15| , 120] for the pure transverse-field Ising model. The method is 
based on the block-spin transformation preserving the high symmetry of the model, and we 
can reproduce the exact critical point and critical exponent v in the pure transverse-field 
Ising chain. This is in contrast with standard real-space renormalization-group approaches 
on the basis of block-spin transformations for quantum systems, which have difficulties in 



quantitatively accurate estimations, see, e.g. 
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FIG. 1. Construction of block spins in one dimension. 



We start by dividing the chain into blocks of two spins as shown in Fig. [TJ The Hamil- 
tonian is also split into intra-block and inter-block parts, 



-J2i-l&2i-l a 2i ~ Ai-lCT. 



2i-l"2i-l' 



(8) 

where spins 2i — l and 2i belong to block i, and spin 2i + 1 belongs to block i + The label 
of block i runs from 1 to N/2. Most importantly, this particular block partition is suited to 
preserve the form of the Hamiltonian under the renormalization-group transformations 
and the key for the success of our calculations. 
The eigenvalues of iy? ntra are degenerate, 
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-(i) 



.(2) 



e {3) = e (4) 



-- ViJ^-i) 2 + (A,-i) 2 . 

The corresponding eigenvectors are 

= a t\ tt) + ar\ |t), = af\ U) + a~\ U), 

|3)i = ar\ ||) - a+\ U), |4)< = o f "| tt) - 4\ It), 

where 



at 



\ 



1 ± 



J: 



2i-l 



(10) 

(11) 

(12) 
(13) 

(14) 



and {| tt), I tt), I tt), I tt)} is the orthonormal basis in the cr^-basis, i.e. a z \ t) = | t), 
o*\l) = -\l). 

We next keep the two lowest lying energy eigenstates |1) and |2), and drop the others, 
|3) and |4), to perform a coarse-graining. This procedure is expected to be effective for the 
study of the ground state. We then replace each block with a single spin representing the 
|1) and \2) states. To this end, we define the projector onto the coarse-grained system as 

N/2 

(15) 



P 



Pi. 



where Pj is the projector, 

P i = (|l)(l| + |2)(2|) i . (16) 

The resulting coarse-grained Hamiltonian is PHP, and in particular the renormalized intra- 
block Hamiltonian is trivially represented with the identity operator lj on block i as 

P i H^ x& P i = e\ 1) l i . (17) 
Terms in the inter-block Hamiltonian are projected as 

Pi(l2i-l®<4 i )Pi = %, (18) 



P i+ i K + i ® l 2l+2 ) P i+1 = =af +1 , (19) 

P (l«-i ® O Pi = t= == =gf, (20) 

where erf is the a-component of the Pauli matrix on block i, or new site i in the coarse- 
grained system. 

The renormalized Hamiltonian is consequently expressed as 

N/2 N/2 N/2 

php = £ 4 ] - E ^ ( 21 ) 

i=l i=l i=l 

with renormalized couplings 

Jj = - J2iJ2i+l =, (22) 

ji _ P2i-\P2i / 2 ox 

4 V(J 2 i-i) 2 + (r 2 i-i) 2 ' 

Note that our transformation preserves the form of the Hamiltonian. In other words, our 
method does not generate additional couplings under renormalization. Other choices of the 
intra- and inter-block Hamiltonians lead to more inconvenient transformations that do not 
preserve the form of the Hamiltonian. 

Let us calculate renormalized A to generate the renormalization group equation, 

i N/2 r 

N l*k * (24) 
= 2A 
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where we have used a property of the periodic boundary condition Jjv/2+1 — A- The 
renormalization group equation A = 2A has a fixed point 



A c = 0, (25) 

and this agrees with the exact transition point [161 ] . Combining the change of A with that of 
the typical correlation length £ typ = £t yP /2 through the scale transformation with the scaling 
factor 2, we obtain the critical exponent 

^typ = 1 (26) 



under Eq. ([7]). This is also equal to the exact value [18]. 

To take untypical effects into account, we have to explore the change of variance of 
(log/ 1 — log J) through renormalization. However, it is difficult to analytically investigate 
it. We therefore study the change of variance by numerical methods. 

In numerical calculations, we first prepare a pool of parameters J and r. Both param- 
eters obey the uniform distributions p(J) = 9(J)9(1 — J) and p(-T) = jr6{r)6{r u — r), 
respectively, where if x > 0, 9(x) = 1, and 9(x) = otherwise. If the upper bound r u of the 
values of transverse fields is equal to 1, the distributions of J and r accord with each other, 
and the system lies on the critical point. The pool contains iV couplings and N transverse 
fields. They are randomly labeled such as Jj and i = 1, 2, . . . , N, and a chain of length 
N is constructed from them. In our calculation, N has been chosen to be 10 6 . 

We next perform the renormalization according to Eqs. ( 122|) and ( 123|) with the periodic 
boundary condition, and the pool is renewed by generating the renormalized couplings and 
fields. Then the size of the system becomes a half of that of the original system. To repeat 
the renormalization on a large system, we add a copy of all couplings and fields in the 
renormalized system to the pool and randomly relabel them. Consequently, the number of 
couplings and fields in the pool is recovered, and the couplings and fields obey the identical 
distribution with that of the pool before the copies are added. It is noted that, since Ji 
and r i+ i in the renormalized system are correlated with each other (Eqs. (I2"2"j) and (|2"5|) ). 
we have to treat Ji and .Tj+i together in the relabeling process, or more specifically, if Ji 
is relabeled as Jj, r i+ % has to be relabeled as -Tj+i- The coupling Ji and field r i+ i in the 
renormalized system are not related with any other parameters, and we thus can rearrange 
the labels as long as the above rule is followed. 

8 
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FIG. 2. Ratio of the square root of the variance of (log-T — log J) after a renormalization to that 
before the renormalization in the random transverse-field Ising model in one dimension. 



Repeating this scheme, we observe the change of the variance of (log-T — log J). To reduce 
statistical errors, we run up to 100 samples. We set the system very close to the critical 
point A c = when the calculation is carried out. More precisely, the upper bound of fields 
r u is set to be 1, and consequently the distributions of couplings and fields are regarded as 
identical ones. Nevertheless these distributions have a bit difference owing to the finiteness 
of the number of couplings and fields, practically. 

A result of the numerical estimate is shown in Fig. [21 where the ratio of the square root 
of the variance of (log T — log J) after a renormalization to that before the renormalization 
is plotted. The ratio is always larger than 1, that is, the distributions of the logarithm of 
parameters keep broadening. This fact demonstrates the existence of an infinite-randomness 
fixed point in the random transverse-field Ising spin chain. Furthermore, the ratio reaches 
a stationary state after several steps of renormalization. Accordingly, the value of critical 
exponent v calculated with Eq. ([6]) also becomes stable, because the ratio of A after a 
renormalization to that before the renormalization is constant, 2. Since randomness is strong 
when the stationary behavior appears, we expect to extract the nature of the system near 
the infinite-randomness fixed point from this stationary behavior. We therefore estimate the 
critical exponent v with the values in the stationary region as 

v = 2.00(7), (27) 

9 



where we have assumed that a renormalized correlation length ought to be a half of the 
original correlation length through the renormalization with scaling factor 2. This value is 
in good agreement with the exact one v = 2 Q . 

Although our real-space renormalization group procedure includes approximations, it 
reproduces the exact critical point, u typ , and v. Our simple scheme correctly reflects the 
physics of the infinite-randomness behavior in particular, which is one of the most peculiar 
features in the random transverse-field Ising spin chain. 



III. TWO DIMENSIONS 

A. Generalization to the two-dimensional model 

We generalize the renormalizat ion-group method to the two-dimensional transverse-field 
Ising model with randomness on the square lattice. This is also a generalization of the previ- 



ous study on the pure model in two dimensions 



15] to the random model. The Hamiltonian 



is 



H 



-£■>. 



_2 _2 

v a i j - 



r £ 



(J- 



(28) 



where spins interact with their nearest neighbors (i, j). The key idea consists of perform- 
ing renormalization-group transformations that preserve the form of the Hamiltonian by a 
projective isometry that preserves the bond algebra (i.e. the algebra realized by the oper- 
ators ofcrj and erf). Exploiting our experience in one dimension, we divide the lattice into 
blocks just as in one dimension (Fig. |3]). Furthermore, we combine the one- dimensional 
block method in horizontal and vertical directions to restore the symmetry of the lattice. 
Specifically, we iterate the renormalization in two directions: first in the horizontal direction 



oo 
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o 
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FIG. 3. Construction of block spins in two dimensions and the two steps of renormalization in the 
horizontal direction and then the vertical direction. 
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and then in the vertical direction (Fig. [3]). 

Now, we redefine the coupling constants for the horizontal direction and the vertical 
direction to distinguish these two quantities in this scheme. The coupling constants between 
spin on and the neighboring spin to the right side is JhU,j) an d that between spin on 
and the neighboring spin to the upper side is J v (i,j), where denotes the location 
of a single site on the two dimensional lattice. 

In the first step of the renormalization (in the horizontal direction) we replace each block 
with a single spin using the same procedure as in the one-dimensional case. We have the 
relation corresponding to Eqs. ( 1TS|) - (I2D|) . 

P(i,j) (l(M-ij) ® 0-( 2iji )) P(ij) = a z (iJ) , (29) 

P (i+hj) ( a hi+l,j) ® 1 (2i+2J)) P(i+l,j) = i = h{2l+1 ' j) ^ ^(i+lj)' ( 30 ) 

y J h(2i+l,j) + r (2i+lj) 

P(i,j) (l(2i-l,i) ® 0-f 2 iJ)) = , = (2t ~ 1,j) = gfij), (31) 

y J h(2i-l,j) + r {2i-X,j) 

where Puj) denotes the projector onto the state space of the block spin, namely the spin 
on in the renormalized system. We find that the z-component of the spin on the right 
spot in a block becomes the ^-component of the block spin, but the ^-component of the 
spin on the left spot in a block becomes the ^-component of the block spin multiplied by 



J hi V 'Jh + P 2 - The renormalized couplings and fields are then written as, 

j Jh(2i,j)Jh(2i+l,j) , . 

M,3) ~ , , r2 =; (32) 

J h(2i+l,j) "T 1 (2i+lJ) 

, J H2i-l,j) Jh(2i-l,j+l) j j ,„o\ 

Jv(i,j) - i 2 ' I 2 2 ^ =J < 2i - 1 >i) + Ji, (2i-i)' \ 66 > 

J h(2i-l,j) + r (2i-lJ) y J h(2i-l,j+l) + F (2i-l,j+l) 

p r(2i-ij)P(2i,j) . s 

= I 2 2 ^ ' 

y J h(2i-i,j) + r (2i~ij) 

In Eq. fl33|) . J v (2ij) is derived from the coupling of two spins on the right spot in each block, 
and the rest is derived from the one on the left spot in the blocks. 

Next the system is renormalized in the vertical direction in the same way as the horizontal 
direction. The coupling constants and the transverse field are now 

J J v(i,2j-1) Jy(i+l,2j-l) J J 

Jh(i,j) — — / = I = = — J h(i,2j-1) + Jh(i,2j), ) 

J v(i,2j-l) + r {i,2j-l) \] J v (i+l,2j-l) + r (i+l,2j-l) 

11 



It is important that our transformations in two dimensions also do not generate extra terms 
as the one- dimensional case, and the form of the Hamiltonian is preserved. In addition, the 
lattice structure is preserved. 

Note that our transformations are not local. After the first step of the renormalization, 
Jh(i,j), Jv(i+i,j), and share Jh(2i+i,j) and r^i+i,j)- Hence, Jh(i,j), Jv(i+i,j), and 

are correlated with each other. Similarly, J/ l (j J+1 ), J„(j +1 j), and -T^+ij+i) share Jh(2i+ij+i) 
and r (2 j + ij + i). Thus, Jh(i,j), Jh(i,j+i), Jv(i+i,j), and are correlated with 

each other. In a column, couplings in the vertical direction, couplings to the left-hand 
spin, and fields are mutually correlated after the horizontal renormalization. The vertical 
renormalization, accordingly, makes couplings in the horizontal direction, couplings to the 
lower-side spin, and fields in a row mutually correlate. After two steps of the procedure, no 
set of couplings and fields are independent from others consequently. 

The renormalization-group transformations are numerically performed as follows. We 
prepare a pool of couplings J and fields r. The pool contains iV couplings in the horizontal 
and vertical directions, respectively, and N fields. In our calculation N is 10 6 . We build a 
two-dimensional partial lattice, which we call a "cluster" hereafter, of couplings and fields 
randomly taken from the pool. The cluster is renormalized according to Eqs. (|32|) - (|37|) . and 
we put the renormalized cluster in a new pool. We set the shape of an original cluster to 
obtain a renormalized square cluster of size L x L. A renormalization process for a cluster 
in the case of L = 2 for example is depicted in Fig. HI The scheme, from building a cluster 
to putting the renormalized cluster in the new pool, for a square cluster of size L x L is 
executed N/ (L x L) times. The new pool is filled by N/ (L x L) renormalized clusters of size 
L x L as a result, and a single renormalization-group transformation for the whole system is 
completed. In the next renormalization step, we construct new clusters of the clusters in the 
pool and repeat the process. In our calculation L is set to be 20. The renormalization for the 
whole system is repeated 15 times, and we run up to 100 samples. The procedure is partially 



based on the calculation method by Nobre 27| in which we randomly take couplings and 



fields from the pool as often as we make a cluster. 
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FIG. 4. Construction of a cluster (left) for the renormalized cluster of size L = 2 (right) in two 
dimensions and two steps of renormalization of the cluster in the horizontal and vertical directions. 



The reason why we have used the clusters in our calculation lies on the nonlocal property 
of our renormalization-group transformations. Although in one dimension we have rear- 
ranged the labels of couplings and fields to repeat the renormalization calculations on the 
large system, we cannot freely do so in two dimensions because the renormalized parameters 
are mutually correlated. It is difficult to fully keep the correlation in the calculation, but 
the cluster procedure takes it into account to a certain extent. We correctly deal with the 
correlation in clusters. Although this procedure still ignores the correlation between the 
clusters on interfaces, the effect of the surfaces is expected to be small if L is large. 



B. Random ferromagnet in two dimensions 

We apply the renormalization-group method to the random ferromagnetic Ising model in 
transverse fields in two dimensions. The couplings J in the pool are uniformly distributed: 
p(J) = 9{J)9{1 — J), and the initial value of the field r is fixed to a constant at all sites. 

We determine the critical point in terms of the parameter (log r — log J) as in the one- 
dimensional case. If the average of (log/ 1 — log J), which is denoted by A, after 15 times of 
renormalization-group transformation is larger than the initial value of A, we conclude the 
system to be in the paramagnetic phase, and otherwise, the system is regarded as in the 
ferromagnetic phase. We obtain the critical value of the transverse field I\ = 0.9115(5) as 



a result. This value is to be compared with another estimate 0.84338(2) [9], which is fairly 
close to our result in consideration of the simplicity of our idea. 

We next observe the change of the variance V of (log/ 1 — log J). The ratio of \/V after 
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FIG. 5. Ratio of the square root of the variance of (log-T — log J) after a renormalization to that 
before the renormalization in the transverse-field random Ising ferromagnet in two dimensions. 
The initial transverse field is set to be 0.920. 



a renormalization to that before the renormalization is plotted in Fig. [H Although the 
result shows a difference form the result in one dimension in that the ratio in the first 
renormalization is smaller than 1, the ratio is always larger than 1 after that. The result 
thus demonstrates the existence of an infinite-randomness fixed point in the two-dimensional 
transverse-field Ising ferromagnet with randomness. This result is same as strong-disorder 
renormalization group jsj and Monte Carlo simulation 



The critical exponent v for the average correlation length is estimated with the relation 

A A.. 2 i/,^_A ; (3g) 



where we have assumed that the renormalized correlation length ought to be a half of the 
original correlation length through the renormalization with scaling factor 2. The symbols A 1 
and V 1 denote A and V renormalized % times, respectively, and A l c means A on the transition 
point after i times of renormalization-group transformations. We regard the average of A 1 
in the cases of the initial transverse field r = 0.911 and 0.912 as A l c . 

The value of A on the transition point depends on the distribution of couplings J. In 
general the distribution changes through the renormalization-group transformations even if 
the system lies just on the transition point. Hence, the value of A on the transition point 
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FIG. 6. Critical exponent v calculated with Eq. (|38p in the transverse- field random Ising ferro- 
magnet in two dimensions. The result in the case F = 0.920, where the system is close to the 
critical point, has a plateau. 



after i times of transformations can be different from the one in the initial distribution. 
Moreover, note that the critical exponent v calculated with Eq. fl38|) can vary with the 
number % of transformations. 

Although the critical exponent v should be estimated near the transition point, insta- 
bilities occur if we try to evaluate it by starting too closely to the transition point due to 
statistical uncertainties. We thus evaluate v with the initial transverse field F = 0.920. The 
results are shown in Fig. [6j For comparison, the result of the case of the initial transverse 
field r = 1.000 is also plotted there. 

We can find a plateau around 4 to 7 for the result of T = 0.920, whereas it does not 
exist in the case of T = 1.000. This plateau is comparable with the stable behavior in the 
one dimensional case. We interpret it as a characteristic property of a critical phenomenon. 
The stable behavior in two dimensions does not continue long owing to the difficulty to keep 
the system close to the critical point due to randomness, whereas we can easily do so in one 
dimension because the fixed point in one dimension is fortunately determined only by the 
average of (log r — log J) and is independent on the distribution of J. This means that in 
one dimension we have to control only the initial value of the average to prevent the system 
from migrating away from the fixed point. Estimating v with values on the plateau (from 4 

15 



to 7 in the horizontal axis), we obtain 



^ = 1.20(6), (39) 

which is consistent with a previous study from a different approach, namely strong-disorder 
renormalization group, v = 1.24(2) [9[. 

C. Spin glass in two dimensions 

We next investigate the spin glass model, where the sign of couplings Jij can take both 
positive and negative values. In our calculations, is independently governed by the 
Gaussian distribution P(Jij) = exp((Jjj — J) 2 /2)/\/27r, where the variance is set to 1. We 
control the average J of the distribution and the uniform transverse field r. 

Let us draw phase boundaries. Which phase the system lies on is determined as follows. 
The paramagnetic phase and the ordered phases, namely the ferromagnetic phase and the 
spin glass phase, are distinguished as in the case of random ferromagnet. Specifically, if 
the average of (log/ 1 — log J) after 15 times of renormalizat ion-group transformations is 
larger than its initial value, the system is regarded as in the paramagnetic phase. Otherwise 
the ferromagnetic or spin glass phase is realized. The criterion for drawing a line between 
the ferromagnetic phase and the spin glass phase is the value of [J] 2 /Vj, where [J] denotes 
the average of Jy and Vj expresses the variance of J^. We conclude that the system is 
in the ferromagnetic phase if [J] 2 /Vj is larger than 1 after 15 times of renormalizat ion- 
group transformations. Otherwise it is in the spin glass phase. We take 30 samples in this 
calculation. 

The resulting phase boundaries are depicted in Fig. [7J The result is not quantitatively 
in precise agreement with a previous study [l^] on the transition point along the line J = 0, 
or more specifically, F c = 1.183(3) in our result and r c = 0.608(4) in 141 ] . It is nevertheless 
important that a definite phase diagram has been obtained, especially with a spin glass 
phase, by the present simple renormalization group with block spin transformations. This 
result implies that our method properly reflects the effect of frustration, which is one of the 
most essential features of spin glasses. 

The ratio of \/V after a renormalization to that before the renormalization is observed 
also in this spin glass model and is plotted in Fig. [HJ which shows the existence of an infinite 
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FIG. 7. Phase diagram of the two-dimensional spin glass in transverse fields. The horizontal axis J 
expresses the average of the Gaussian distribution of couplings Jij , and the vertical axis expresses 
the magnitude of the transverse field r. The symbols, P, SG, and F, denote the paramagnetic, 
spin glass, and ferromagnetic phases, respectively. 

randomness fixed point in the present system. This result supports a conjecture on the 
possible existence of an infinite-randomness fixed point in spin glasses inferred inderectly 
(i.e. without explicit calculations) from the strong-disorder renormalization group approach 
to the random ferromagnet system [5(. 

We next calculate the critical exponent v for the average correlation length. We focus 
ourselves on the transition point of J = 0. As in the case of random ferromagnet, the 
estimation is executed slightly away from the actual transition point F = 1.183(3) to avoid 
instabilities. Specifically, the initial transverse field is set to be r = 1.195. The exponent v 
calculated in the same way as in the random ferromagnet is shown in Fig. [91 To compare 
the critical region with another case away from criticality, results with the initial transverse 
field r = 1.300 is also plotted there. We can find a difference between the two results that 
the value of v for T = 1.195 reaches a plateau after four times of transformations, and the 
value increases again after ten times of transformations. Since this behavior corresponds to 
the plateau in the random ferromagnet, we accordingly evaluate u with values in this region 
(from 4 to 9 in the horizontal axis), 



v = 1.21(9). 



(40) 
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FIG. 8. Ratio of the square root of variance of (log r — log J) after a renormalization to that before 
the renormalization in the two-dimensional spin glass in the transverse field r = 1.195. 
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FIG. 9. Critical exponent v calculated with Eq. 1138 j) in the two-dimensional spin glass in transverse 
fields. The result in the case of F = 1.195, where the system is close to the critical point, has a 
plateau. 



For comparison, if we use the values in the same region as the case of the random ferromagnet 
(from 4 to 7 in the horizontal axis), we have v = 1.19(8). 

In addition to having an infinite-randomness fixed point, the value of v estimated with 
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FIG. 10. Three steps of renormalization in three dimensions in the regular order A and the reverse 
order B. 



our method for the spin glass is in good agreement with that for the random ferromag- 
net (Eq. (13"9"j) ). This fact possibly suggests that the random ferromagnetic Ising model in 
transverse fields and spin glass in transverse fields would belong to the same universality 
class. 

IV. THREE DIMENSIONS 

A. Generalization to the three-dimensional model 

We go on to the three-dimensional case. We generalize our renormalization-group scheme 
to three dimensions. The transformations are iterated in the horizontal, vertical, and third 
directions consecutively (Fig. [10]). We define the coupling constants between spin (i,j,k) 
and the neighboring spin along the third direction (i,j,k + 1) as Jt(i,j,k), and similarly for 
Jh(i,j,k) an d J v (i,j,k) for the horizontal and vertical directions. 

In the first step of renormalization (in the horizontal direction) the parameters change as 




Jh(2i+l,j,k)Jh(2i+2,j,k) 



(41) 





(43) 



(42) 
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Note that the coupling constants of the vertical and third directions are changed under the 
same rule. Carrying out analogous transformations in the vertical direction and then the 
third direction after this first step, we obtain the parameters of the system renormalized in 
the three directions. 

This scheme, however, yields anomalous results due to the imbalanced treatment between 
three directions. As the transformations for three directions are iterated, the magnitude of 
coupling constants of the third direction J t rapidly becomes small in comparison with those 
of the other two directions despite the fact that the original system has no anisotropy. The 
scheme thus needs some corrections. 

To resolve the anisotropy problem, we further renormalize the system in the reverse order, 
along the third, vertical, and then horizontal directions after the regular order described 
above (Fig. [10]) . The six steps, the regular order and then the reverse order, is regarded as 
a single transformation of scaling factor 4. This modified procedure succeeds in rendering 
virtually isotropic results. This is the same process introduced in the previous study for the 
non-random system jl5| . 




Numerical calculations are implemented in the same way as in the two-dimensional case. 
We generate a pool containing N couplings in the horizontal, vertical, and third directions, 
respectively, and N fields, where N is 10 6 . The cluster method is used also in three dimen- 
sions with the cluster of size L x L x L, and L is set to be 5. We repeat the (six-step) 
renormalization-group transformations 5 times, and 100 samples are run. 

B. Random ferromagnet in three dimensions 

We first examine the random ferromagnetic Ising model in transverse fields on the cubic 
lattice. The initial distribution of couplings J in the pool is p(J) = 9 (J) 9(1 — J) identical 
to the two dimensional case, and we control the magnitude of the uniform field r at the 
initial condition. 

We estimate the transition point by comparison of the initial value of A and the final 
value of A after 5 times of transformations. The result is T = 1.266(2). 

The ratio of the square root of the variance V of (log r — log J) after a renormalization 
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FIG. 11. Ratio of the square root of variance of (log-T — log J) after a renormalization to that 
before the renormalization in the transverse-field random Ising ferromagnet in three dimensions. 
The initial transverse field is set to be 1.269. 

to that before the renormalization is plotted in Fig. [TTJ The stability of the ratio through 
renormalization is not obvious, but at least the result demonstrates the existence of an 
infinite-randomness fixed point. 

The critical exponent v for the average correlation length is also computed with the 
relation 



where we regard the average of A 1 in the cases of the initial transverse field T = 1.266 and 
1.267 as A l c . We can find a plateau in the result of r = 1.269 most clearly near the critical 
point (Fig. [T2l . Evaluating v with the values on the plateau (from 2 to 4 in the horizontal 
axis), we have 



A i +l _ M +l 



A 1 - A\ 



(45) 



4 




v = 0.92(4) 



(46) 



which is close to 0.97(5) of the corresponding result by strong-disorder renormalization 
group [lOj, 111 ]. 
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FIG. 12. Critical exponent v calculated with Eq. (|45p in the transverse- field random Ising ferro- 
magnet in three dimensions. 

C. Spin glass in three dimensions 

We apply our renormalization-group scheme to the spin glass in transverse fields on 
the cubic lattice. As in the case on the square lattice, we use the Gaussian distribution 
P(Jij) = exp((Jjj — J) 2 /2)/\/27r, and we control the average J of the distribution and the 
uniform transverse field r. 

We first draw phase boundaries under the same rule as in two dimensions, provided that 
phases are determined after five transformations. The resulting boundaries are depicted in 
Fig. [T31 It is significant that our simple renormalization-group method verifies the existence 
of the spin glass phase. 

Next, critical properties are investigated. We specifically treat the critical point for J = 0. 
The existence of an infinite-randomness fixed point is confirmed by the observation of the 
change of \/V through the renormalization-group transformations (Fig. [T4"l) . The critical 
exponent v is also calculated. We use the average of A 1 in the cases of the initial transverse 
field r = 1.613 and 1.612 as A l c . A plateau is clearly recognized in the result of F = 1.619 
(Fig. IT5|) . and we determine v with the values on the plateau (from 2 to 4 in the horizontal 
axis) , 

v = 0.94(3). (47) 
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FIG. 13. Phase diagram of the three-dimensional spin glass in transverse fields. The horizontal 
axis J expresses the average of the Gaussian distribution of couplings Jij in the model, and the 
vertical axis expresses the magnitude of the transverse field r. The symbols, P, SG, and F, denote 
the paramagnetic, spin glass, and ferromagnetic phases, respectively. 

This value is close to that of the random ferromagnet model in three dimensions (Eq. fj46|) ). 



V. CONCLUSION 

We have studied the transverse-field Ising model with randomness in finite dimensions 



by the real-space renormalization group method introduced in previous studies [15l . 120 ] . 
Our method reproduces exact results for the transition point and critical exponent v in 
one dimension in spite of the existence of randomness. Moreover, our generalization of the 
method to higher dimensions has been shown to be effective not only in the pure model [ijj] 
but also in the random model. In fact, we have obtained the values of the critical exponent 
v in the two- and three-dimensional random ferromagnetic Ising model close to those from 
the strong-disorder renormalization group approach [9l4ll|. 

Most remarkable results in our study concern two- and three-dimensional spin glasses in 
transverse fields. First, we could draw phase diagrams in the J — T plane, which has not been 
done so far to the best of our knowledge. Second, we have concluded the existence of infinite- 
randomness fixed points from indefinite amplifications of the distribution of parameters 
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FIG. 14. Ratio of the square root of variance of (log-T — log J) after a renormalization to that 
before the renormalization in the three-dimensional spin glass in transverse fields with r = 1.619. 
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FIG. 15. Critical exponent v calculated with Eq. (|45p in the three-dimensional spin glass in 
transverse fields with F = 1.619. 



in logarithmic scale. This observation supports the conjecture inferred from the case of 
random ferromagnet under strong-disorder renormalization group 5] but is in conflict with 
a relatively old numerical study 14|. Furthermore, the estimated exponent v in spin glasses 
in transverse fields is very close to that of the corresponding random ferromagnet. Thus it 
is naturally expected that these models belong to the same universality class, which is also 
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a highly non-trivial result . 

The validity of our method is not readily obvious, because it drops higher-energy eigen- 
states in the block Hamiltonians. We nevertheless believe that consistencies of our results 
with previous studies, wherever applicable, justify our procedures to extract the essential 
features of random quantum systems. It would be necessary to further improve the method 
to achieve quantitatively better results. 
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